Signatures of hadron-quark mixed phase in gravitational waves 
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We calculate stellar oscillations including the hadron-quark mixed phase considering the finite 
size effects. We find that it is possible to distinguish whether the density discontinuity exists or not 
in the stars, even if one will observe the gravitational waves of the fundamental mode. Additionally, 
the normalized eigenfrequencies of pressure modes depend strongly on the stellar mass and on the 
adopted equation of state. Especially, in spite of the fact that the radius of the neutron star with 
IAMq , which is standard mass, is almost independent from the equation of state with quark matter, 
the frequencies of pressure modes depend on the adopted equation of state. Thus, via observing the 
many kinds of gravitational waves, it will be possible to make a restriction on the equation of state. 
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PACS numbers: 04.40.Dg, 97.60.Jd, 26.60.Kp 



I. INTRODUCTION 



In order to detect the gravitational waves, which are oscillations of spacetime itself, several ground-based detectors, 
such as LIGO, VIRGO, TAMA300, and GEO600, are in operation and there are some projects to build the next 
I , generation of detector In addition to the ground-based detectors, it is also considering to launch detectors in 
Q ■ space like LISA @ and DECIGO 3] . Since the permeability of gravitational waves could be extremely strong, one 
can expect to see raw information of the wave sources via gravitational waves. On the other hand, the most promising 
sources of gravitational waves might be supernovae and mergers of binary compact objects, i.e., the gravitational 
waves are related to the compact stars, around which the gravitational field should be strong. So, the direct detection 
of gravitational waves could enable us not only to collect the astronomical data and to reveal the true properties of 
dense matter (34Tl| , but also to prove the gravitational theory in the strong-field regime [H, [TH . 

In fact, an attempt to estimate the stellar parameters, such as mass, radius, and equation of state (EOS), via 
their oscillation properties is not a brand-new idea. In astronomy, the helioseismology has been already established 
and one could see the interior information of our sun through its oscillation properties. Since the late 1990s, it has 
been suggested possible to reveal the compact star properties with observing the oscillation spectra, which is called 
gravitational wave asteroseismology [4,5]. Furthermore, the detailed analysis of the emitted gravitational waves might 
permit us to determine the radius of accretion disk around supermassive black hole [14j or to see the magnetic effect 
' during the stellar collapse [Hj]. 

With respect to the neutron stars, the density in the vicinity of the stellar center could become much more than 
the standard nuclear density, which is around po rs 0.17 fm -3 . Since such high density cannot be realized on the 
earth, the detailed matter properties in the neutron stars are still unknown. However, this means conversely that 
the neutron stars can be good candidates to know the matter properties in the extreme high density region, where 
it is suggested that the exotic components of matter, such as hyperons, meson condensates, and quark matter, could 



appear [T(j. The existence of these exotic matter changes the EOS and neutron star structures dramatically |17H19| . 
Namely, as mentioned the above, observing the stellar oscillations and/or the corresponding gravitational waves will 
tell us the information about the matter properties of neutron stars. Especially, it might be impossible to probe the 
true properties of matter deep inside the star by any other experiments. 

The hyperons are considered to appear at around 2 — 3po, if the nuclear matter would be in beta-equilibrium [20l. l2lj] . 
On the other hand, there are still many uncertainties with respect to the hadron-quark phase transition, e.g., the EOS 
of quark matter or a deconfinement mechanism. The presence of quark matter inside the compact objects mi ght play 
an important role in the astronomical phenomena, for example the backbending effect from phase transition |22| . l23j , 
connection to gamma-ray bursts [24| , the gravitational wave bursts [2pl , the gravitational radiatio n l26r]29i , the energy 
release during the collapse from neutron stars to quark stars [3(1 l3l| . the neutrino luminosities |32l434| . and cooling 
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with quark matter [35l439l |. Although the simplest model with the quark matter would be the Maxwell construction, 
around critical density to appear the quark matter the mixed phase could exist, where baryon number and electric 
charge should be conserved. Generally, the properties of the mixed phase depend strongly on the electromagnetic 
interaction and the surface tension, which effects are called "thefinite-size effects" . Due to the finite-size effects, the 
mixed phase is composed of the nonuniform pasta structure (l8l.[l9|. However, it has not been clear how to distinguish 
the finite size effects by observing the astronomical phenomena. 

In this article, we study the gravitational waves emitted from compact stars with the hadron-quark mixed phase 
considering the finite size effects. Up to now, there exist a few studies from the gravitational wave asteroseismological 
point of view, which focus onprobing the density discontinuity in the high density region with using specific frequencies 
of gravitational waves (e.g. , @, [|[ ) . However, no one investigates the effects of the mixed phase between the hadron and 
quark matter on the specific frequencies. So, in this article, preparing the neutron star models with the mixed phase 
or with density discontinuity between the hadron and quark matter phases, we will calculate the eigenfrequencies 
associated with the gravitational waves, where as a first step we adopt the Cowling approximation, i.e., the metric 
perturbation will be neglected. Then, as varying the stellar properties systematically, we will see the dependence of 
the existence of quark matter on the oscillation frequencies considering finite size effects. 

This article is organized as follows. In the next section, we describe the equation system to construct the neutron 
star models and the adopted EOS in this article, and show some stellar models concretely. In section Hill we derive 
the perturbation equations with Cowling approximation. With appropriate boundary conditions, the problem to solve 
becomes the eigenvalue problem. Then the obtained oscillation spectra will be shown in section HV1 At last, we make 
a conclusion in section [V] In this article, we adopt the unit of c = G = 1, where c and G denote the speed of light 
and the gravitational constant, respectively, and the metric signature is (— , +, +, +). 



II. NEUTRON STAR MODELS 



The equilibrium configurations of non-rotating relativistic stars are spherically symmetric solutions of the well- 
known Tolman-Oppcnheimer-Volkoff (TOV) equations. The metric can be expressed as 

ds 2 = -e 2 *dt 2 + e 2A dr 2 + r 2 d9 2 + r 2 sin 2 dd(f> 2 , (2.1) 

where $ and A are metric functions with respect to r. A mass function m(r) are defined as m(r) = r(I — e~ 2A )/2, 
which satisfies 

dm , ,„ „, 

— = 47rr 2 p, (2.2) 
dr 

where p is the energy density, while the TOV equations to determine the distributions of the pressure P(r) and metric 
function $(r) are described as 

dP , ,d<S> , . 

<P + P)-T, (2-3) 



dr dr 
d$ to + 47rr 3 P 



dr r(r — 2m) 



(2.4) 



To close the equation system, one needs an additional equation, i.e., the equation of state (EOS). 

In this article, we adopt the EOS with the hadron-quark mixed phase with hyperons considering finite size effects 
according to [f8], [lj| . Our EOS for hadron is in the framework of the nonrelativistic Brueckner-Hartree-Fock a ppr oach 
including hyperons such as £~ and A 40]. It is not clear now whether the £~-iV interaction is repulsive or not [4l[|42j|. 
while we use here a weak but attractive interaction. It would be interesting to see how our results are changed by 
using the other E~-A interactions, and we will discuss it in the future work. For the comparisons, we also adopt the 
EOS composed of only nucleon. We called them as "hyperon EOS" and "nucleon EOS" in this article. 

For the quark phase, we adopt the MIT bag model. It should be noticed that the adopted EOS in this article is not 
simple MIT bag model but more sophisticated models suggested in the previous articles [431 EH ■ Assuming massless 
u and d quarks and s quarks with the current mass of m s = 150 MeV, we set the bag constant B to be 100 MeV 
fm- 3 . 

For the mixed phase, we assume non-uniform structures, so-called "pastas" . In practice, the structures such as 
droplet, rod, slab, tube, and bubble are considered. We use the local density approximation for particles using the 
Wigner-Seitz cell. In order to construct such pasta phase, we put a sharp boundary with a constant surface tension 
parameter between the hadron and quark phases. Then the Gibbs conditions are imposed, where one needs to solve 
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the chemical equilibrium among particles in two phases consistent with the Coulomb potential, and a pressure balance 
consistent with the surface tension in self-consistently. Although the knowledge of the value of the surface tension at 
the hadron-quark interface is very poor, some theoretical estimations have done and they suggest that the value of 
surface tension is around a ~ 10 — 100 MeV/fm 2 [45|, HfjJ. Since one can see that the models with <r> 40 MeV/fm 2 
are almost same as that with a = 40 MeV/fm 2 , in this article we consider only a range of 10 < a < 40 MeV/fm 2 . As 
the extreme case, we also consider the EOS with the Maxwell construction, which has the sharp discontinuity of the 
density between 5.93 x 10 14 and 8.82 x 10 14 g/cm 3 . Note that this discontinuity appears at strong surface tensions 
considering the finite size effects, e.g., a > 70 MeV/fm 2 in our previous study [la]. Moreover, we take into account 
another extreme case of the EOS with the bulk Gibbs condition, which appears at the zero surface tension limit [l8| . 
Finally, in order to determine which structure is most favored in the mixed phase, we compare the energy among the 
pasta structures. It notes that we do not take into account the anti-particles and muons in this article, because their 
effects should not be so important. 

Then the above EOSs should be connected to the hadronic EOS proposed by Negele and Vautherin 47J when the 
density become lower than around 10 14 g/cm 3 . In Table U the components of each EOS adopted in this article are 
summarized, while Fig. [T] shows the relationship between the energy density p and the pressure P in the higher 
density region for the adopted EOSs in this article. Solving the TOV equations with such EOSs, we can get the stellar 
properties as Fig. [21 It should be noticed that the maximum masses of the stellar models composed of EOS including 
the quark matter are almost independent from the value of a, which are around 1.4Mq. At the end of this section, 
we should mention about the maximum mass of neutron stars. Recently two pulsar mass measurements appeared, 
which are well above 1.4M©, i.e., M = 1.97M Q for PSR J1614-2230 48] and M = 1.667M© for PSR J1903+0327 [H- 
Since it is impossible to explain these evidences with the EOS including the quark matter adopted in this article, we 
should derive more realistic EOS and examine the stellar oscillations as a future work. 



TABLE I: Components of each EOS adopted in this article and the references. 



EOS 


components 


nucleoli 


nucleon 


hyperon 


nucleon, hyperon 


bulk Gibbs 


nucleon, hyperon, quark 


with pasta phase 


nucleon, hyperon, quark 


Maxwell 


nucleon, hyperon, quark 




FIG. 1: Relationship between the total energy density including masses (p) and the pressure (P) for the EOSs adopted in this 
article. 



III. PERTURBATION EQUATIONS 

In order to determine the specific oscillation frequencies in the neutron stars, in this section we present the per- 
turbation equations for nonradial oscillations of spherically symmetric neutron stars. Especially, in this article we 
consider only stellar oscillations with Cowling approximation, where the fluid would oscillate on a hxed background 
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FIG. 2: Mass-radius relation for neutron stars with several EOSs. 



metric. Namely, the spacetime will be frozen such that the metric perturbation should be neglected (5g^ v — 0). Thus, 
with Cowling approximation, our study limits to the modes related to the fluid perturbations, i.e., /, p, and g-modes, 
while we cannot see the emission of gravitational waves associated with the so-called iu-modes which correspond to 
the spacetime oscillations. It should be emphasized that even with this simple approximation one can see qualitatively 
the features for oscillation frequencies of emitted gravitational waves. 
The fluid Lagrangian displacement vector is given by 

C = (e~ A W, -Vde, -V sin" 2 68$) r-' 2 Y tm , (3.1) 

where W and V are functions with respect to t and r while Yp m is the spherical harmonic function. Then the 
perturbations of the four- velocity, Su^, can be written as 

8u» = (Q,e- K d t W,-d t Vde,-d t V sin' 2 6d^)r- 2 e-*Y tm . (3.2) 

With these variables, the perturbation equations describing the fluid oscillations can be obtained by taking a variation 
of the energy-momentum conservation law, i.e., S (\7 V T^) = 0. This equation reduces to V V 5T^ U = with Cowling 
approximation. The explicit forms with p — r,6 are 



+ ^ ( 1 + ^) [e- A W + £{£ + 1)V] - ?-±J^&e- K W = 0, (3.3) 



r 2 \ dp 

(p + P)e- 2 *V + 2f [e- A W + £{£ + l)V] + ^e~ A W = 0, (3.4) 

where primes and dots on the variables denote the partial derivative with respect to r and t, respectively. 7 is the 
adiabatic constant defined as 

_ fd\nP\ _ nAP 

1 ={~dh^) s ~PA^' (3 - 5) 

where n is the baryon number density and A denotes the Lagrangian variation. The Lagrangian variation of the 
baryon number density, An, is determined by the relation as An/n = — v[, 3 '£ fe — Sg/g, where Vj^ and g are the 
covariant derivative in a three dimension with metric g^ v and the determinant of g^ u , respectively. Since in this 
article we adopt the Cowling approximation, the second term is neglected. Then the Lagrangian variation of n can 
be described as 

^ = -e"^ - ^V. (3.6) 

Assuming a harmonic dependence on time, the perturbative variables will be written as W(t, r) — W(r)e and 
V(t, r) = V(r)e lu)t . Additionally, calculating the combination of the form d[Eq. (|3.4|) ]/dr — [Eq. (|3.3p ] and substituting 
Eq. (13.41) again, one can get the simple equation as 

V = 2&V -e A ^. (3.7) 
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Thus, from Eqs. (|3 .4[) and (|3.7I) . one can obtain the following simple equation system for the fluid perturbations; 



W'=^ [wVe A - 2 * V + &W] - £{l + l)e A V, (3.8) 



V = 2<P'V-e A %. (3.9) 



In order to solve this equation system, we have to impose appropriate boundary conditions at the stellar center 
(r = 0) and at the stellar surface (r = R). With these boundary conditions, the problem to solve becomes an eigenvalue 
problem for the parameter uj. From the above equation system, one can find the behavior of W and V in the vicinity 
of stellar center as W(r) = Cr e+1 + 0(r e+3 ) and V(r) — —Cr e /£ + 0(r i+2 ), where C is an arbitrary constant. On the 
other hand, the boundary condition at the stellar surface is the vanishing of the Lagrangian perturbation of pressure, 
i.e., AP = 0. Since AP could be expressed that AP — -fPAn/n from Eq. (|3.5|) . with the help of Eqs. (|3.6I) and 
(13.81). the condition of AP = becomes as 



2 r 2 e A - 2 *V + &W = 0. (3.10) 



Furthermore, if one would consider the stellar models with density discontinuity, one has to prepare the additional 
junction condition at the surface of discontinuity, which are the continuous condition for W and AP Q. These 
junction conditions can be rewritten with the variables W and V as 

W+ = W-, (3.11) 



p2 $ 

V + = 



u 2 R 2 



( [u 2 R 2 e~ 2 *V- + e- A &W-} - e - A &W+ ) , (3.12) 

I P+ + P J 



where R g denotes the position of the density discontinuity, and W-, V_, and p_ are the vales of W, V, and p at 
r = R g — while W + , V + , and p + are the values of W, V, and p at r = R g + 0, respectively. 



IV. OSCILLATION SPECTRA 



In this section we examine the stellar oscillations on the stellar models shown in Sec. HU Especially, we focus on the 
stellar models whose mass is in the range of 0.5Mq < M < M max and at O.IMq intervals, where M max is maximum 
mass for each EOS. Namely, the masses of the stellar models we adopt in this article are 0.5 < M/Mq < 1.3 for hyperon 
EOS, 0.5 < M/Mq < 1.8 for nucleon EOS, and 0.5 < M/M Q < 1.4 for the other EOSs. As mentioned the above, 
the stellar models with Maxwell EOS have the density discontinuity, if the central density is larger than 8.816 x 10 14 
g/cm 3 . That is, for Maxwell EOS, the stellar models with 0.7 < M/Mq < 1.4 have the density discontinuity, while 
those with M/Mq = 0.5 and 0.6 do not have the density discontinuity and such stellar models are same as those with 
hyperon EOS (see Fig. [2]). 

When neutron stars oscillate, many kinds of gravitational waves are radiated. If the stars are spherically symmetric 
and without density discontinuity inside the star, which might be the simplest model, the fundamental (/), pressure 
(p), and spacetime (w) modes arc excited, where / and p modes are gravitational waves related to the fluid oscillations 
while w modes correspond to the oscillations of spacetime itself. If the stars are spherically symmetric and with density 
discontinuity, the additional oscillation modes, i.e., the g modes, are excited as well as /, p, and w modes. The g 
modes are also gravitational waves associated with the fluid oscillations. In this article, we will see qualitatively the 
gravitational waves related to the fluid oscillations because the Cowling approximation is adopted in our analysis. 
Thus, as shown in Fig. [3J the stellar models with the adopted EOS except for Maxwell EOS have / and p mode, 
while those with Maxwell EOS, whose masses are more than 0.7M©, have /, p, and g modes. 

Before discussing the / and p modes, we pay attention to the g mode for the stellar models with Maxwell EOS. As 
noted before, g mode is excited due to the existence of density discontinuity. Therefore, one could know the existence 
of density discontinuity inside the neutron stars, if the g mode gravitational waves will be observed. Actually, since 
the typical frequency of g mode is in the range from a few hundred Hz up to kHz, such gravitational waves could be 
observed by using the ground-based gravitational waves detectors. From Fig. [31 one can observe that the frequency 
of g mode is almost independent from the stellar mass, which is around 1.73 Hz. However, we find that g mode 
frequency can be expressed well as a function of stellar compactness M/R (see Fig. @|, such as 

wM = 0.3130 (—] +0.0103. (4.1) 
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FIG. 3: The first a few eigenfrequencies for the stellar models with (a) nucleon and (b) Maxwell EOSs are plotted as a function 
of the stellar mass M/Mq, where the frequency / is defined as / = lu/2h. As mentioned in text, for the case of the stellar 
models with the other EOSs, the kinds of excited eigenfrequencies are same as to the case of the stellar model with nucleon 
EOS, i.e., they have / and pi modes, where i = 1, 2, 3, ■ • ■ . On the other hand, for the case of the stellar model with Maxwell 
EOS, which has 1st order phase transition in the density, the additional eigenfrequency are excited, which is g mode. 
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FIG. 4: The normalized eigenvalue ui of g mode for the stellar models with Maxwell EOS are plotted as a function of the stellar 
compactness M/R. 



Practically, this empirical formula can expect the g mode frequency with less than 0.6% accuracy. So, via observing 
the g mode gravitational waves, one could know the stellar properties. Or, with the observation of redshift parameter, 
which is connected to the stellar compactness directly, one could make a restriction on the stellar mass. 

On the other hand, it is well-known that the frequency of / mode can be connected to the stellar average den- 

1 /2 

sity [M/R 3 ) , which could be physically explained by considering the relation between the sound speed and the 
propagation time of the fluid perturbation inside the star. In fact, Andersson and Kokkotas derived the empirical 
formula for the / mode frequency by calculating in the stellar models with the several realistic EOSs, which did not 
include the EOS considering the quark matter, and found that the obtained / mode frequencies in their article are 
subject to this empirical formula almost independently of the adopted EOS [5]. While, in Fig. [5J we show the / mode 
frequencies for the stellar models with adopted EOSs in this article. At a glance, one can observe that the behavior 
of / mode frequencies for the stellar models with Maxwell EOS is quite different from the others. These deviations of 
the results for Maxwell EOS become up to 36% from those for nucleon EOS and 27% from those for the other EOSs. 
This means that one could possible to know the existence of the density discontinuity even by observing the / mode 
gravitational waves as well as g mode ones. Additionally, Fig. [5] shows that the / modes for the stellar models with 
pasta EOS are similar to those for the stellar models with the other EOS without density discontinuity. 

Regarding to the p modes, we plot the frequencies of p\ and pi modes as functions of average density in Fig. [S] 
From these figures, one can see that pi modes frequencies are almost independent from the EOS even if that includes 
the density discontinuity. However, one can see the dependence of p mode frequency on the EOS if we make the 

1 /2 

figure of the normalized eigenvalues with average density, uj [R? /M) , as a function of the stellar mass (see Fig. 
[7]). The most interested point in this figure is that the normalized eigenvalues for the stellar models with IAMq 
depend strongly on the adopted EOS, in spite of the fact that the stellar shapes are almost independent from the 
adopted EOSs including quark matter, i.e., R — 9.42 — 9.67 km. Moreover, in this figure, the dependence of the 
normalized eigenvalues of p\ mode on EOS is different from that of pi modes. Thus, with the help of the observations 
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FIG. 5: With several EOSs, the frequencies of / modes are plotted as a function of the stellar average density (Mj R 3 ) 1 ^ 2 , 
where // is defined as // = loj /2tt. The left panel corresponds to the results with nucleon, hyperon, bulk Gibbs, and Maxwell 
EOSs, while the right panel focus on the pasta EOSs. 



of stellar mass, it could be possible to distinguish the EOS by observing the several kinds of oscillation modes. At 
last , it should be also noticed that the normalized eigenvalues of / modes for the stellar model with Maxwell EOS 
are obviously different from the other stellar models as well as Fig. [5] 




(M/R 3 ) m {MIR 3 ) 




FIG. 7: The normalized eigenvalues of / (left panel), pi (middle panel), and pi (right panel) modes are plotted as a function 
of the stellar mass M/Mq, where the marks in figures are corresponding to those in Fig. [5] 
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V. CONCLUSION 

We study how to distinguish the finite effects on the hadron-quark mixed phase by observing the gravitational 
waves, for which we derive the perturbation equations of neutron stars and obtain their eigenfrequencies. 

We find that one could know the existence of density discontinuity inside the star via observing the gravitational 
waves of not only g mode but also / mode. Note that this discontinuity comes from the instability of the mixed 
phase due to the strong surface tension. Additionally, it is possible to see the stellar properties by observing the 
g mode frequency, since such frequency can be expressed well as a function of the stellar compactness. If the EOS 
do not include the density discontinuity, it might be difficult to distinguish the EOS only by observing the / mode 
frequencies. However, the normalized eigenfrequencies of p modes depend strongly on the EOS even if the EOS do 
not including the density discontinuity, although the raw frequencies of p mode is almost independent. Thus, with 
the help of the observation of stellar properties, it could be possible to make a restriction on the stellar EOS. 

In this article, for simplicity, we assume the Cowling approximation, which restricts our examination to only stellar 
oscillations. This means that we should do a more detailed study including the metric perturbations. Via this 
type of oscillations, we could obtain the additional information, such as the damping rate of gravitational waves, 
and combining with results shown in this article would provide more accurate constraints on the stellar properties 
and/or the stellar EOS. Furthermore the stellar magnetic field might play an important role. For example, the quasi- 
periodic oscillations have observed in the giant flare and these phenomena are believed to be oscillations of strong 
magnetized neutron stars Considering the effects of stellar magnetic fields, it might be possible to obtain the 

further information. 

Additionally, although we focus only on neutron star matter (T = MeV and Y Ve — 0) in this article, in order 
to study the proto neutron stars, we should take into account other effects, e.g., the thermal effects on the stellar 
oscillations [5l[, and the effects of temperature and/or neutrino trapping on the pasta structures [l9l l52l [53j. 
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